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Specific lieat and resistivity measurements were performed on polycrystalline samples of the solid- 
solution Ya:Lui-a;Ni2B2C in order to determine thermodynamic properties such as the specific- 
heat difference AC, the thermodynamic critical field Hc(T), as well as the upper critical field 
Hc2{T). These properties were analyzed within the Eliashberg theory including anisotropy effects, 
yielding electron-phonon coupling anisotropy parameters (o^) ranging between 0.02 and 0.03 for the 
whole series, and Fermi velocity anisotropy parameters of (b^) = 0.245-0.3. Excellent agreement 
between theory and experiment was achieved for these parameters, the Sommerfeld constant 7 and 
model phonon spectra determined from specific heat measurements. An analysis of the previously 
investigated boronitride La3Ni2B2N3_i for comparison revealed the electron-phonon anisotropy to 
be of great significance in describing its thermodynamic properties and the calculations yielded 
(al) ~ 0.08 and (bl) ~ 0.245. The Tc-behavior within the series Ya;Lui_a;Ni2B2C is discussed in 
terms of he density of states at the Fermi level A*'(0). 



I. INTRODUCTION 

The transition-metal borocarbide superconductors 
i?Ni2B2C with transition temperatures comparable to 
those of the A-15 compounds (e.g. i? = Lu with Tc ~ 
16.5 K, NbsGe with Tc ~ 23 K) are a subject of broad 
interest for research on intermetallic superconductors. 
Siegrist et aln reported the crystal structure of the 
i?Ni2B2C superconductors to be a filled version of the 
ThCr2Si2-type structure stabilized by the incorporation 
of carbon, where Ni2B2 layers are separated by RC lay- 
ers. In the related compound La3Ni2B2N3_5, three LaN- 
planes separate the Ni2B2 layers. Despite of the lay- 
ered structure reminiscent of the cuprate superconduc- 
tors, the electronic structure of the single i?C-layer boro- 
car bides as well ^J^i-fhe triple LaN-layer boronitride is 
three-dimensional.H^B Nickel-site substitutions on both 
the borocarbides and the boronitride revealed similar 
electronic-, properties of the bands related to the 3d- 
electrons.cl 

At a first glance, the thermodynamic properties of 
La3Ni2B2N3_5 seem to be close to the weak coupling 
BCS-predictionSjB but standard single-band BCS-theory 
cannot explain the pronounced upward curvature of the 
upper critical field Hc2{T) close to Tc- If applied to 
the borocarbides, it fails to describe both the thermo- 
dynamic properties and the uraper critical field obtained 
from experiment. Shulga et aln analyzed the upper crit- 
ical field Hc2{T) of LuNi2B2C and YNi2B2C in terms of 
the Eliashberg-theory using an isotropic two-band model. 
An isotropic single band model cannot reproduce the pos- 
itive curvature near apparently, because of the domi- 
nant role of anisotropy effects in this system. Recently, 



Dugdalc et alU presented an experimental investigation 
on the Fermi surface of LuNi2B2C, revealing that it con- 
sists of three sheets. m 

Freudenberger et alu studied the solid solution 
Ya;Lui_2;Ni2B2C and showed, that Tc exhibits a mini- 
mum at about x ~ 0.5 with a Tc-depression of about 
IK. They pointed out, that this feature cannot be de- 
scribed alone by disorder effects with the residual resis- 
tivity ratio as a measure, and tentatively supposed the 
electron phonon coupling strength A to be the origin of 
this behavior. 

We present in this paper investigations of the ther- 
modynamic properties and of the upper critical field 
Hc2{T) of the series Ya;Lui_a:Ni2B2C and the boroni- 
tride La3Ni2B2N3_5. The critical temperature, the spe- 
cific heat, and the temperature dependence of the up- 
per critical magnetic field were measured for all samples 
of the series Ya;Lui_a;Ni2B2C and the results, includ- 
ing previous measurements of La3Ni2B2N3_5 (Ref. ||), 
were analyzed using an anisotropic model of the s-wave 
Eliashberg-formalism. Considering the change in mass 
from Y to Lu in the series and different lattice properties 
of both systems, we used model phonon spectra calcu- 
lated from the phonon contribution to the specific heat 
and found these to be sufficient in describing both the 
thermodynamic properties and the upper critical field 
of La3Ni2B2N3_5 and Y3:Lui_a;Ni2B2C. An analysis of 
the coupling strength and electronic density of states in 
Yj,Lui_2;Ni2B2C gives an insight on the origin of the dip 
in Tc(a;). Sec. II presents details of sample preparation 
and measuring techniques and Sec. Ill shows the theoret- 
ical background used for the analysis. In Sec. IV, an anal- 
ysis of the experimental data in terms of the anisotropic 
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Eliashberg theory is discussed and possibilities for the 
Tc-reduction are presented. Conclusions are drawn in 
Sec.V. 



II. EXPERIMENTAL 

Polycrystalline samples of Yj;Lui_2:Ni2B2C were syn- 
thesized on a water-cooled copper grove by high fre- 
quency induction melting under argon atmosphere. The 
starting materials are rare earth ingots (R = Y, Lu) 
(Strem chemicals, USA: 99.9%), Ni ingots (Strem chem- 
icals, USA: 99.9%), crystalline boron (Starck, Germany: 
99.5%) and carbon ingots (Starck, Germany: 99.99%). 
A good homogeneity was obtained by performing the fol- 
lowing three steps of sample preparation: (i) Nickel and 
boron were melted two times, (ii) Rare-earth and car- 
bon ingots were alloyed together while compensating the 
carbon loss (1-2%) during the remelting stage, (iii) The 
two precursor alloys NiB and RC were melted together. 
The buttons were broken and remelted for at least 12 
times, and finally annealed at 1020°C for at least one 
week in evacuated quartz tubes. The phase purity of 
the samples was checked at room temperature by apply- 
ing X-ray diffraction (XRD) in a Guinier-Huber camera 
using germanium as an internal standard. The X-ray 
photographs of the samples, indexed on the basis of the 
tetragonal crystal structure of LuNi2B2C (space group 
lA/mrnm), show that the samples of the solid solution 
Y2;Lui_,2.Ni2B2C are almost single phase with traces of 
Yi_2;Lu2;B2C2 aud Yi_a;Lua;Ni4B. According to the x- 
ray line intensities and optical micrograph investigations 
the total amount of these secondary phases is about 2-5% 
in all samples investigated. 

Specific heat measurements in the temperature range 
1.5K-160K and magnetic fields up to 9 T were carried 
out on 1-2 g samples employing a quasi-adiabatic step- 
heating technique. The sample holder consists of a thin 
sapphire disc (m ^ 0.2 g) with a strain gauge heater and 
a CERNOX temperature sensor. The field calibration 
of the latter has been performed against two GaAlAs 
resistivity thermometers and a capacitive SrTiOs sensor. 
Four point resistivity measurements were performed in 
fields up to 9 T in order to determine Hc2 (T) by means 
of a midpoint criterion. The transition widths according 
to a 10% and 90% criterion are indicated by error bars 
in the figures containing the upper critical field obtained 
from the experiment. 

The characterization of the polycrystalline 
Ya;Lui_2;Ni2B2C samples with respect to their residual 
resistance po and the room temperature residual resis- 
tance ratio (RRR) turned out to be rather problematic, 
because our well annealed samples are rather brittle and 
partly textured. Thus, we had difficulties to obtain re- 
producible and reliable values of the RRR ranging be- 
tween 7 and 41. We obtained a very high RRR value of 
41 for LuNi2B2C (twice as large as typical single crystal 



values, see e.g. Ref. ||) but a much lower value of just 12 
for YNi2B2C although its thermodynamic mean transi- 
tion temperature Tc = 15.56 K (see below) is even slightly 
higher jthan Tc = 15.4 K of a high quality YNi2B2C single 
crystaJia with RRR= 37. The upper critical fields Hc2{T) 
match to each other despite of the different RRR values 
(compare Fig. 11 of this workjjsath Fig. 4 of Ref. |l^). In 
the latter paper Nohara et allEi reported that the Som- 
merfeld value 7 of clean limit YNi2B2C and LuNi2B2C 
exhibits a ^/H dependence which changes to a linear re- 
lation 7 cx -ff in the dirty limit. We note, that we found 
an approximate 7 cx \/H dependence in all our samples 
including La3Ni2B2N3_5, although the RRR shown in 
Ref. I is only 3, similar to the "dirty" Y(Nio.8Pto.2)2B2C 
for which 7 oc iJ was reported. It is obvious that the 
grain boundaries of the polycrystalline samples strongly 
reduce the RRR value especially in La3Ni2B2N3_5. Fur- 
ther arguments for the extrinsic origin of the large resid- 
ual resistance po of La3Ni2B2N3__5 will be given in con- 
text with analysis of the experimental data in terms of 
the anisotropic Eliashberg theory (see Sec. VIE). 



III. THEORETICAL BACKGROUND 

Anisotropy effects have to be considered in the descrip- 
tion of the experimental results for Hc2{T) in order to 
obtain a satisfactory agreement pbetween theory and ex- 
periment. Daams and Carbottetll applied the separable 
model introduced by Markovitz and KadanoffE3 to de- 
scribe an anisotropic electron-phonon interaction spec- 
tral function within Eliashberg theory: 



[a'-F('^)]k,k' = (1 + ai,)a^F{Lu){l + at') : 



(1) 



where k and k' are the incoming and outgoing quasi- 
particle momentum vectors in the electron-phonon scat- 
tering process and at is an anisotropy function describ- 
ing the deviation of the anisotropic spectral function, 
[a^F(ijj)]k k', from the isotropic one, a^f (w), in the di- 
rection of k. It has the important feature that its Fermi 
surface average (ak) — and as anisotropy effects are 
rather small it is sufficient to keep only the mean-square 
anisotropy (a^) as the important anisotropy parameter. 
The theory of Hc2{T) for anisotropic polycrystalline 

,me was de- 
It employs 



superconductors in a separable model sc. 
veloped by Prohammer and Schachinger 
the separable ansatz for the anisotropy of the electron- 
phonon interaction and the ansatz 



(1 + MM, 



(2) 



which describes the anisotropy of the Fermi velocity; 6k is 
an anisotropy function defined the same way as ak. The 
upper critical field is then described by the set of equa- 
tions (27-30 in Ref. |l^ and thermodynamic properties 
of the borocarbides were calculated using Eqs. (33-37) of 
Ref. A^-band models have been extensively studied 
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as a tool to describe|-|anisotropic features of supercon- 
ductors. Sbulga et alu applied the two-band model Hc2 
equationfJlj to their experimental results and were able 
to describe the experimental data with the corresponding 
equations. Thus, the separable model employed in this 
work can be described in its simplest form by a Fermi 
surface split into two half-spheres of equal weight 

P{a) = 5{-a)/2 + d{a)/2, (3) 

with radii r ± a, if r js the radius of the equivalent 
isotropic Fermi sphercllj Using the Fermi surface har^ 
monies (FSH) notation introduced by Allen,t£l DaamsEZi 
observed that this separable model was equivalently de- 
scribed by a restriction to zeroth-order FSH in each of 
the two subregions of the Fermi surface. According to 
her work, the separable model applied in our analysis 
corresponds to a two-band model in which the two Fermi 
surface regions have equal weight. It is of course also 
possible to define different weights for the two regioiia, 
thus changing the distribution function (||) , but DaamstZl 
demonstrated that in the case of weak anisotropy the in- 
fluence of different weights in a separable model is of neg- 
ligible significance for the thermodynamics of anisotropic 
superconductors. 

In the case of the upper critical field any deviation from 
the equal weight configuration causes i?c2 (T) to approach 
the isotropic case as is demonstrated in Fig. |l|. The de- 
viation function in the insert of Fig. ^ 

Dh.. (T/T,) = Dh^, (t) = §^ - 1 (4) 

demonstrates the deviation of the upper critical field of 
an anisotropic system Hc2,a{T) from the upper critical 
field of an isotropic, equivalent system Hc2,i(T). The nu- 
merical result for equal weights (1:1) was fitted to the 
upper critical field of LuNi2B2C to fix the parameters 
(a^), (6^), and (vf)- The mean Fermi velocity (vf) and 
its anisotropy parameter (6^) were used to fit the exper- 
imental data near Tc, and (a^) was changed to describe 
Hc2{T) at lower temperatures. With these parameters, 
we calculated the upper critical field for different weights 
of the Fermi sheets (l:^) to investigate the change in the 
behavior of Hc2- Obviously, Hc2.a approaches Hc2.i as 
one of the two sheets becomes dominant in weight and 
this is indicated by a flattening of Dh^2 i^) ■ 

The relative signs of at and 6k in the same Fermi- 
surface sheet is also of importance to the analysis of Hc2 ■ 
With same signs of Ok and fek in same sheets, Hc2.a{T) 
is reduced compared to Hc2,i{T) (Fig. ^). Opposite signs 
within same sheets give rise to an enhancement of Hc2,a 
over Hc2,i- In this case, the reduction of Hc2,a{T) due 
to the influence of (&^) is compensated by the electron- 
phonon interaction anisotropy. Bandstructure calcula- 
tions can give a good estimate for the weight and the rel- 
ative signs of Ok and 6k in the corresponding regions of 
the Fermi surface. Dugdale et aln were able to show that 
the Fermi surface of LuNi2B2C consists of three sheets 



without giving the weights belonging to the Fermi surface 
regions. The third one is a small electron pocket centered 
at r, which compared to the other two sheets can be re- 
garded to have a negligible effect on Hc2 if considered in 
the theoretical calculations. 

With the knowledge of the Fermi surface consisting of 
two dominant sheets and the shape of the experimental 
upper critical field favouring opposite signs of ak and 6k 
within the same sheets, we used equal weights for the 
two Fermi surface regions and opposite signs of Ok and 
6k in the following analysis. Numerical results for Hc2{T) 
showed that in the case of both anisotropy functions hav- 
ing the same signs, the mean Fermi velocity {vf) needed 
to fit the experimental data deviates significant^ from 
the value determined from the plasma frequencyll3. The 
same effect on {vf) is observed with different weights of 
the Fermi surface regions. |— ■ 

Finally, impurities are treated in Born's limito which 
assumes the impurities to be randomly distributed and 
to be of dilute concentration. In such a limit impurities 
are characterized by a scattering rate i"*" which is pro- 
portional to the impurities' concentration. Their main 
effect is the smearing out of the electron-phonon interac- 
tion anisotropy resulting in a slight reduction of TcJi3 an 
enhancement of Hc2 at low temperatures, and a reduc- 
tion of the high temperature upward curvature of Hc2 as 
was demonstrated by Weber et a/.o for Nb. 

IV. RESULTS & DISCUSSION 

A. Results of the specific heat measurements 

The low-temperature specific heat of Yj;Lui_a;Ni2B2C 
with x = 0, 0.1, 0.5, 0.8, and 1 is presented in Figs. 

An external field of 9 T was applied to suppress su- 
perconductivity in order to determine the normal state 
heat capacity C„ = Ce + Cph = ■yT + f}T^, where 7 
is the Sommerfeld parameter and P is related to the 
low-temperature value of the Debye temperature. From 
the specific heat data in the normal and superconduct- 
ing state (Figs. ^ and ^ we obtained the specific heat 
difference AC — Cg — Cn and used it to calculate 
further quantities as the thermodynamic critical field 
Hc{T) and its deviation function D{T) (Eqs. ^^), 
which are compared to numerical results obtained with 
Eliashberg theory as discussed below. The extrapola- 
tion of the normal state specific heat to T ^ shown 
in Fig. ^ with the range of extrapolation starting at 
^ 30 revealed the Sommerfeld constants of the series 
Y^Lui_:rNi2B2C ranging from 7 = 18.5(2) mJ/molK^ 
for X = 0.8 to 20.6(2) mJ/molK^ for x = 0. The 
values for the other compounds within the series are 
19.7(3) mJ/molK^ {x = 0.1), 18.7(3) mJ/molK^ [x = 
0.5), and 19.7(2) mJ/molK^ {x — 1). A comparison of 
the Sommerfeld constants of the border phases [x = 0,1) 
with those of Ref. |2^ reveals about ^ 1 mJ / mol 
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higher values of 7 in this work, with about the same 
for LuNi2B2C but a higher one for YNi2B2C compared 
to previous results. The difference in these properties 
may be explained by the homogeneity range achieved 
by differeiit methods used to prepare polycrystalline 
materials. EJ 



B. Spectral function and coupling strength 

The large mass difference between Y and Lu causes a 
significant change of the low energy phonon frequencies 
in Y2;Lui_j;Ni2B2C while high energy phonon modes are 
hardly affected by the rare earth substitution. The lat- 
ter was shown by inelastic neutron spectroscopy (INS^ 
results on Y- and LuNi2B2C reported by Gompf et a/.E3 
Similar conclusions were also obtained from the lattice 
heat capacity of i?Ni2B2C [R = La, Lu and Y) analyzed 
by means of a modeL phonon density of states (PDOS) 
in a previous paper .E2I The low energy part of our model 
spectra were well confirmed by INS while the high en- 
ergy part is, of course, just a crude description which 
does not give the real details. Nevertheless, the mo- 
ments of these spectra are in remarkable agreement with 
those calculated from the INS results. As the following 
discussion of superconducting properties in terms of the 
Eliashberg formalism requires primarily details of the low 
energy PDOS we evaluated the variation of the low en- 
ergy phonon modes within the series Y2;Lui_2;Ni2B2C by 
a similar procedure as explained in Ref. i.e. by fitting 
the normal state specific heat measured from 2 to 160 K 
with a model PDOS, F{uj), which is of the form 

F{uj) = m eiu;)e{nu - lj)^ + 

9{Lj~LUo)9inE3-^))- (5) 

In this formula 9{x) is the well-known step function; 
luq — 4.3 meV in order to avoid a finite spectral density 
at zero energy; ai are the widths of the Gaussians and 
gi are the weights of these contributions (i.e. the num- 
ber of the contributing phonon branches). The Debye 
spectrum (oc ui^) represents the spectral weight of the 
three acoustic branches and the Gaussian contributions 
account for the 15 optical branches. The free param- 
eters are the Debye cut-off and the Gaussian peak 
positions fl Ei , respectively. In order to fix a high energy 
limit of the model spectrum, the Gaussian contribution 
with the highest energy (ri^,) was cut off at the peak 
position and 53 is accordingly doubled in equation (||). 
The parameters applied in Eq. (|^) are summarized in 
Table |. The only difference with respect to the model 
spectra shown in Ref. is the splitting up of Ei with 
gi = 1.5 into two parts ^Eo ^nd U,Ei with go — 0.3 and 



gi = 1.2, respectively, which was motivated by the shape 
of F{w) obtained by INS on YNiaBaC (see Ref. ||). On 
the other hand we now fixed fJ^;., to 103.4 meV (equiva- 
lent to 1300KxfcB) in order to keep the number of pa- 
rameters small. Only for Yo.5Luo.5Ni2B2C we applied a 
slightly modified spectrum with g^ = 0.7 and gi = 0.8 
compared to 0.3 and 1.2 given in Table |, because of the 
fifty to fifty ratio of Y and Lu in this sample. 

Adjusting the model spectra to the lattice heat capac- 
ity via 

r°° ( 
C,h{T)^R du;F{u;) . ^j'. , (6) 
Jo smh / 

we obtained the model spectra shown in Fig. |^. This pro- 
cedure provides a satisfactory description of the Cp{T) 
data from 2K up to 160 K (up to 300 K for the border 
phases) . 

The phonon contribution to the specific heat yields in- 
formation on F{uj) rather than on a'^F{u!). Therefore, 
one is led to introduce an assumption on the form of the 
electron-phonon coupling function a^ioj). Junod et a/.Ej 
suggested a function (w) ^ oj" in & first approximation 
and found the exponent to be s = —1/2 for a large num- 
ber of A15 compounds by comparing the experimentally 
determined moments of a'^F{ijj) with tunneling data. We 
calculated the thermodynamic properties of the borocar- 
bides using s — Q, —1/2, and —1, and found the exponent 
s = —1/2 to be the best choice for the theoretical analysis 
as in the case of the A15 compounds. It was impossible to 
describe the experimental data using s = because the 
theoretical deviation function D{t) already lies below the 
the experiment for an isotropic system, as a consequence 
of the large gaussian contribution above the low temper- 
ature Debye frequency ^ and is expected to get more 
negative for an anisotropic system. In the case of s = — 1 
the calculated deviation function is significantly larger 
than the experimental one and a consistent description 
of the thermodynamic properties and Hc2 is impossible 
because of the large anisotropy parameter (a^) necessary 
to describe the thermodynamic properties. 

The spectra were cut off at the local minimum of the 
model phonon spectra at w ~ 68meV. Starting at this 
energy, the phonon spectra determined by Gompf et alM 
using inelastic neutron-scattering techniques have a wide 
interval, where the PDOS is zero. The next optical con- 
tributions to F{uj) start at ^ 100 meV, and it is not very 
likely that electrons still considerably couple to phonons 
at these energies. 

With the coupling function having the form a^^uj) = 
auj~^^^, the model spectra shown in Fig. ^were weighed 
with a function uj^^^'^ and were rescaled with a constant a 
to give the corresponding critical temperatures of the 
samples for a fixed pseudopotential of /i* = 0.13. The 
model spectral functions give electron-phonon interaction 
mass enhancement factors 

X^2r,.^ (7) 

Jo w 
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between 1.02 {x = 1) and 1.22 {x = 0). A does not behave 
as the critical temperature and descends monotonically 
from LuNi2B2C to YNi2B2C as a function of x, without 
showing a minimum, which indicates that it cannot ac- 
count alone, if at all, for the depression of Tc within the 
series. 



C. Thermodynamic properties and Hc2{T) 

With the set of parameters given in Table ||, the ther- 
modynamic properties were calculated for all samples. 
The number of atoms in a volume of 1 cm^, n^, was cal- 
culated from lattice parameters a and c of the tetragonal 
structure. The Sommerfeld constants obtained from the 
experiment, ^exp, are slightly greater than the theoreti- 
cal values, 7c, needed to fit the experimental results for 
the thermodynamic critical field Hc{T). The latter is 
calculated from the experimental data by integrating the 
difference between the zero-field and 9T measurement, 
i.e. any normal state contribution to the zero-field heat 
capacity cancels out in AC = - C„ = C(0 T) - C(9 T) 
which is used in the analysis below. Nevertheless, a nor- 
mal state contribution, which is eliminated in AC, will 
strongly influence the experimental result for the tem- 
perature dependence of the electronic specific heat in the 
superconducting state Ces(T), obtained by subtracting 
the phonon contribution from the zero-field data. We 
note, that the 7-value due to a small amount of impu- 
rity phase obtained by extrapolating the zero field data 
in a C/T vs. plot to T ^ yields e.g. for YNi2B2C 
Jimp — 0.4(2) mJ/molK^ which is significantly smaller 
than jexp — Ic ^ l.lmJ/molK^. For all other samples, 
Jimp is found to be even closer to zero, but jexp—Jc is still 
about 1 mJ/molK^ (see Fig. |^). Thus, the average differ- 
ence of about 5% may be attributed to an intrinsic effect, 
namely to the contribution of the small electron pocket, 
for which de Haas-vaji Alphen measurements performed 
by Terashima et a/.CJ revealed the superconducting gap 
to be much smaller than on other parts of the Fermi 
surface. For a further discussion of the temperature de- 
pendence of Ces(r) see below (Sec. IV F). 

In the description of the thermodynamics and the up- 
per critical field of Yj:Lui_a;Ni2B2C the model spectra of 
Fig. H were used in the numerical solution of the Eliash- 
berg equations in the clean limit (impurity scattering rate 
=0). This is of course only a rough approximation 
and the corresponding scattering rate t'^ should be de- 
termined from a Tc reduction with increasing concentrjEtj 
tions of Born limit momentum scattering impurities.cj 
As it was not possible to perform this type of experiment 
with the present samples we investigated the influence 
of impurity concentration theoretically and realized that 
in order to stay within the constraints of the known ex- 
perimental data had to be rather small and therefore 
we regarded the clean limit to be an acceptable approxi- 
mation. As the coherence length is rather short in these 



materials the clean limit is always a good approximation. 
The results of our theoretical study of different spectral 
functions a^i^(cj) and of the influence of impurities on 
the thermodynamics as well as on Hc2 will be the topic 
of a separate, forthcoming publication. 

Fig. ^ depicts the comparison between experimen- 
tal results and theoretical predictions for LuNi2B2C. 
The top frame shows the specific heat difference be- 
tween superconducting and normal state AC(T) — 
Cs{T) ~ C„(r). Theoretical calculations were done 
for an isotropic ((a^) ~ 0) and an anisotropic case 
((^k) ~ 0.03), and results with anisotropy parameters 
between these values lie in between the numerical results 
of Fig. ^(a-c). Considering impurity scattering would 
bring the results of an anisotropic calculation closer to 
the isotropic case. 

The thermodynamic critical field Hc{T) is calculated 
from the free energy difference /S.F 

Hc{T) = V^AF/^io , (8) 

and the specific heat difference AC is related to AF 
through 

Strong coupling and anisotropy effects can be observed 
in the deviation function of the thermodynamic critical 
field 

Error bars of the experimental deviation function result 
from experimental error bars of the specific heat data and 
primarily from the extrapolation of Hc{T) to T — > 0. 
Anisotropy effects push the deviation function towards 
the weak coupling regime as shown in Fig. ^(c). Even 
with the large error bar attached to the experimental 
data, it becomes obvious that the BCS result is far off 
with its minimum at about —3.7%. 

The upper critical field Hc2{T) was fit to the ex- 
perimental data by describing the positive curvature at 
Tc with the Fermi velocity anisotropy parameter (&^) 
at a fixed mean Fermi velocity (wf), and by changing 
the anisotropy parameter (a^) of the electron phonon 
coupling strength to describe the behavior of the ex- 
perimental curve at lower temperatures. The value of 
{vf) was taken from the experimental result for the 
plasma frequency htOpi ~ ^47re2w|,A^(0)/3 = 4.0 eV (see 
Ref. [1^ ) which results in a mean Eermi velocity {vp) of 
~ 0.28 X 10^ m/s for LuNi2B2C.Q A reasonable fit of 
the experimental data with this value is achieved with 
(a^) = 0.02 and (6^) = 0.25. According to the thermo- 
dynamic properties, it would as well be possible to apply 
a higher anisotropy parameter, i.e. (a^) — 0.03. Error 
bars of the plasma frequency were not given in Ref. [l^ , 
and a change of 0.01 x 10^ m/s in [vp) is likely to be 
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within the experimental error range. For an anisotropy 



parameter of 



0.03, one has to reduce (6^) and 



enhance {vp) to compensate the change as depicted in 
the bottom frame of Fig. ||. The error bars of the exper- 
imental upper critical field result from the width of the 
transition in the resistivity data. 

Figs. ^, |l^, |ll|, and |l^ show the results for the sam- 
ples with X — 0.1, 0.5, 0.8, and 1. Thermodynamic 
properties and Hc2 of the whole series can be described 
with anisotropy parameters (a^) between 0.02 and 0.03, 
which is relatively small compared to the Fermi velocity 
anisotropy (5^) ranging from 0.245 to 0.3 and is com- 
parable to the electron-phonon anisotropy parameters of 
LaAl2 (Ref. || 



) = 0.0335) 

and is smaller than the anisotropy parameter of NbsSn 



= 0.01), Nb (Ref. 14 



(Ref. 



0.08). The Fermi velocity anisotropy pa- 



rameters of Ya;Lui_j;Ni2B2C are obviously higher than 
those of LaAls ((&|) = 0.16), Nb {{bl) = 0.118), and 
NbsSn {{bD — 0.13). In comparing the upper critical 
field of YNizBaC to that of Yo.sLuo.sNiaBaC it is obvi- 
ous that the positive curvature of Hc2 at Tc is rather pro- 
nounced in both samples. If the reduction had been 
caused by impurity scattering the positive curvature of 
Hc2 at Tc would have been gradually reduced with in- 
creasing impucity concentration as was demonstrated by 
Weber et alM in their analysis of polycrystalline Nb. 
Such an effect can certainly not be observed in our case. 



D. Tc-reduction in the series Ya;Lui_i,Ni2B2C 

The reduction of Tc within the series cannot be de- 
scribed with impurity scattering in terms of Eliashberg 
theory because this would require randomly distributed 
scattering centers. A substitution of Lu with Y or vice 
versa produces inhomogeneities of the periodic potential 
only on the rare earth lattice site, so that a random dis- 
tribution can only be expected with small substitutions 
up to a: = 0.1 . If the Tc-reduction should be explained 
by inelastic impurity scattering expressed by the scat- 
tering rate , it is restricted to the condition that the 
electronic properties of the system, as e.g. A^(0), remain 
unchanged, which is not the case in Yj;Lui_2,Ni2B2C (See 
Fig. p^ ). Thus, the Tc-rcduction within the series cannot 
be attributed to inelastic impurity scattering expressed 
hyt+. 

The density of states N{0) (Fig. ^3|) was calculated us- 
ing the bare Sommerfeld constant 7* = 7/(1 + A) and 
the number of atoms per cm'^, n^, in the approxima- 
tion N(0) cx UA^* ■ Unlike the electron-phonon coupling 
factor A, 7 behaves like Tc as a function of x, therefore 
accounting for a similar behavior of N{0). Hence, as A 
decreases with decreasing Lu-content, iV(0) shows a min- 
imum between x = and 1. Obviously, the behavior of Tc 
results from an interplay of the coupling strength of elec- 
trons to phonons and the density of states at the Fermi 
level. The latter is most likely due to a band broadening 



arising from the i?-ion size mismatch which washes out 
the maximum of the electronic density of states at the 
Fermi energy. The N{0) values of the unalloyed com- 
pounds YNi2B2C (4.14 states/eVf.u.) and LuNi2B2C 
(3.94 states/eVf.u.) agree closely with the values of 4.31 
and 4.06 states/eVf.u., respectively, repopted from den- 
sity functional calculations by Divis et alEi 



E. Thermodynamics and Hc2 of La3Ni2B2N3_5 

The model spectrum of La3Ni2B2N3_5 (Fig. |l^) was 
obtained in a similar manner as the spectra of the boro- 
carbides (Ref. |2^) and is cut off at the local minimum at 
~ 50 meV. After applying a coupling function (w) ~ 
with s = —1/2, the spectrum was rescaled with a con- 
stant factor to give the critical temperature Tc — 11.73 K 
and a fixed pseudopotential n* = 0.13. This results in 
a coupling strength of A = 1.02. With the Sommerfeld 
constant 7 — 24mJ/molK^, the thermodynamic proper- 
ties were calculated for an isotropic and an anisotropic 
clean limit case with (a^) = 0.08 (Fig. We note, 

that the clean limit approximation may, of course, be 
questioned because of the rather high residual resistance 
repoted in Ref. || where all data used for the present anal- 
ysis have been measured on one individual sample. Sup- 
posing the high po — 9 /if2cm to be intrinsic, we had to 
consider an impurity scattering rate i"*" ~ lOmeV which, 
however, did not allow to find any reasonable set of pa- 
rameters describing the experimental results on the up- 
per critical field and specific heat. In the view of the 
present Eliashberg approach, these data are compatible 
with scattering rates < 1 meV indicating that the high 
po — 9 /ificm of La3Ni2B2N3_5 comes from scattering on 
the grain boundaries rather than from scattering inside 
the grains which were shown by high resolution electron 
microscopy to be practically defect frecB 

The upper critical magnetic field Hc2(T) shown in 
Fig. lU was fitted with the anisotropy parameter ob- 
tained from the analysis of the thermodynamic proper- 
ties, (al) = 0.08, a mean Fermi velocity (vp) = 0.225 x 
10^ m/s and its anisotropy parameter of = 0.245 . A 
fit to the experimental data with an anisotropy parame- 
ter of (a^) — 0.06 is included for comparison. The mean 
Fermi velocity obtained in the analysis lies between the 
root mean square anisotropics = 0.292 x lO^m/s 

and -v/^jl) = 0.148 x lO^m/s calculated by Singh and 
Pickett .0 In contrast to the borocarbides a consistent de- 
scription of the thermodynamic properties and the up- 
per critical field can also be achieved with s = 0, which 
results in an electron-phonon coupling factor A = 0.87, 
anisotropy parameters (a^) = 0.06 and (6^) = 0.27, and 
a mean Fermi velocity of (vp) = 0.2 x lO^m/s. Differ- 
ent coupling factors A are obtained with different cou- 



pling functions a {w) 



(s = 0,-1/2) because, af- 



ter rescaling the spectral functions to obtain the crit- 
ical temperature of the system and a pseudopotential 
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/i* = 0.13, the low-energy contributions to a^F(u;) are 
enhanced with s = — 1/2 compared to the case of s = 0. 
The coupUng factor A obtained with the couphng func- 
tion 0^(10) ~ Lo^ is in good agreement with the estimate 
A = 7/7f)s — 1 — 0.85. In general, both descriptions 
(s = 0,-1/2) require a higher value for the anisotropy 
parameter (a^) compared to the borocarbides in order to 
describe the experimental data, because the agreement 
between theory and experiment improves with increasing 
(a^) as is clearly demonstrated by the comparison with 
the specific-heat difference, the thermodynamic critical 
field, and the deviation function of La3Ni2B2N3_5 shown 
in Fig. |l^. The rather large anisotropy of (a^) ~ 0.08 fur- 
ther explains the temperature dependence of the upper 
critical field and iJc2(0) of ^ 8T comparable to that of 
LuNi2B2C, since (a^) has an enhancing effect on Hc2- 

F. Comparison of Ya;Lui_i;Ni2B2C and La3Ni2B2N3_i 

The most distinctive feature of the borocarbide and 
boroniti'id|e superconductors pointed out in previous 
papersBca is the temperature dependence of the elec- 
tronic specific heat in the superconducting state: While 
Y- and LuNi2B2C exhibit an almost cubic temperature 
dependence of Ces{T), La3Ni2B2N3_5 shows an approx- 
imately exponential behavior of Ces{T) close to the BCS 
predictions. Thus, La3Ni2B2N3_5 was supposed to be a 
weak-coupling BCS superconductor, because other char- 
acteristics of the superconducting state of the boronitride 
like the thermodynamic ratios (e.g. AC/jTc ~ 1.4) and 
the deviation function D{t) (see Fig. |l^) were also found 
to be much closer to the BCS predictions than those of 
the borocarbides. The above analysis of the thermody- 
namic properties and the upper critical field, however, 
demonstrates that the difference between the thermody- 
namic ratios and the deviation function of the borocar- 
bides and the boronitride arises from a larger anisotropy 
of the electron-phonon interaction in the latter rather 
than from a smaller coupling strength (A ~ 0.87 or 
1.02 for La3Ni2B2N3_5 compared to A ~ 1.02 - 1.22 
for Ya;Lui_a;Ni2B2C). The different temperature depen- 
dence of the heat capacity in the superconducting state 
- cubic in case of the borocarbides - cannot be explained 
by our Eliashberg calculations which fit the specific heat 
difference AC = Cs — C„, but not the absolute Ces{T), 
of e.g. YNi2B2C depicted in Fig. ^ as a semilogarith- 
mic plot of Ces/lTc versus the inverse reduced temper- 
ature Tc/T. As the gap obtained by the above calcu- 
lations is free of nodes, it yields an exponential tem- 
perature dependence for Ces{T) which becomes linear 
in the semilogarithmic representation. This is in con- 
trast to the experimental result for YNi2B2C which ex- 
hibits a pronounced curvature in this plot indicative for 
a power law. In section IV C we noted, that a satis- 
factory description of the thermodynamic properties of 
the borocarbides and in particular their thermodynamic 



critical field is only achieved when the calculations are 
performed with gamma values of about 5% smaller than 
the experimental values. By adding a linear contribu- 
tion corresponding to the 5% of the normal state gamma 
value (1.1 mJ/molK^) omitted before, an improvement 
of the temperature dependence of the calculated Ces (T) 
is achieved (dashed line in Fig. |l^). These additional 
electronic contributions may be attributed to a normal 
state impurity phase and/or to contributions froin the 
small electron pocket, for which Terashima et alt3 pro- 
posed a significantly reduced gap. The latter argument 
is supported by the specific heat results-obtained on a 
single crystal YNi2B2C by Nohara et alrB which are in- 
cluded in Fig. |l^ for comparison. For temperatures Tc/T 
between 1 and 5 (i.e. T = 15.5-3 K) both experimental 
data sets are in rather good agreement, but clearly de- 
viating from our calculations. Below about 3K an ex- 
ponential temperature depeadence can be resolved from 
the data of Nohara et alH. This may be attributed 
to the small gap on the electron pocket which can be 
identified as the-,third Fermi surface sheet detected by 
Dugdale et alB't3 The third Fermi surface sheet was not 
considered in our calculations because its contribution 
appeared to be negligible in — C„. The exponential 
low temperature behavior of Ces in particular for the 
polycrystalline YNi2B2C is covered by a linear electronic 
contribution from small amounts of a secondary phase of 
about 0.4(2) mJ/molK^. Although this contribution is 
significantly reduced in the samples of the pseudoquater- 
nary system, the overall feature of the semilogarithmic 
plot is preserved. Accordingly the conclusions concern- 
ing the reduced gap on the third Fermi surface sheet can 
also be drawn for Ya;Lui_a;Ni2B2C. 

V. CONCLUSIONS 

Model phonon spectra obtained from specific heat 
measurements, rescaled to give the critical tempera- 
tures of the system and for a Coulomb-pseudopotential 
value of /i* = 0.13, resulted in coupling strengths A 
ranging from 1.02 for YNi2B2C to 1.22 for LuNi2B2C. 
We were able to show that a consistent description of 
the thermodynamic properties and the upper critical 
field of the series Y2;Lui_a;Ni2B2C and the boronitride 
La3Ni2B2N3_5 within the Eliashberg theory is achieved 
only, if anisotropy effects of the electron-phonon coupling 
and of the Fermi velocity are included. Excellent agree- 
ment between theory and experiment was achieved in the 
analysis of the series Ya;Lui._a;Ni2B2C with anisotropy 
parameters (a^) = 0.02-0.03 and = 0.245-0.3. The 
model spectral functions for La3Ni2B2N3_^ yielded cou- 
pling factors A = 1.02 (s = -1/2) and A = 0.87 (s = 0) 
comparable to those of Ya;Lui_j;Ni2B2C. The thermody- 
namic properties and the upper critical field are well de- 
scribed with (al) = 0.08 and (bl) = 0.245 (s = -1/2) or 
(al) = 0.06 and (bl) = 0.27 (s = 0). For La3Ni2B2N3_5 
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it is of importance to include anisotropy in order to de- 
scribe its upper critical field and - in contrast to the boro- 
carbides - its thermodynamic properties. 

Coupling and impurity effects are not capable 
of explaining the Tc-behavior within the series 
Ya;Lui_j;Ni2B2C. The behavior of 7 as a function of x 
suggests N{0), aside from A, to be the reason for the 
minimum in Tc{x), which is confirmed by an approxi- 
mate calculation of iV(0) from 7. The approximately 
cubic temperature dependence of Ces(T) of the borocar- 
bides cannot be accounted for by the current Eliashberg 
calculations due to a small gap of the third Fermi sur- 
face sheet. This contribution, however, is suppressed in 
AC and therefore the Eliashberg calculations are in good 
agreement with the related experimental thermodynamic 
properties. 
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TABLE I. Parameters of the model spectrum 



cont. 










Q.E3 (cut-off) 




free 


free 


free 


free 


103.4 meV (fixed) 


ai (meV) 




1.6 


1.6 


12.1 


25.9 


gi 


3 


0.3 


1.2 


8.5 


5 
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TABLE II. Parameters used for the calculation of thermo- 
dynamic properties of Yi,Lui-a;Ni2B2C. The mass enhance- 
ment factor or clcctron-phonon coupling factor, A, the number 
of atoms in a certain volume, ua, and the Sommerfeld con- 
stant used in the theoretical description, 7c, are compared 
within the series Ya;Lui-xNi2B2C. 



X 


Tc [K] 


A 


riA [10^^ /cm^] 


7c [mJ/mol K^] 





16.01 


1.220 


9.405 


19.7 


0.1 


15.56 


1.182 


9.382 


18.9 


0.5 


14.90 


1.058 


9.263 


18.1 


0.8 


14.68 


1.018 


9.220 


18.0 


1 


15.45 


1.017 


9.175 


18.6 



FIG. 1. Influence of different weights for the two 
Fermi-surface sheets on the upper critical field Hc2{T). 

The numerical result in the case of equal weights (1:1) 
is a fit to the upper critical field of LuNi2B2C with 
(o^) = 0.02, (b^) = 0.25, and (vf) = 0.280x10^ m/s. Dn^^it) 
is depicted in the small insert. 



FIG. 2. Influence of difi^erent signs of Ok and fek on the 
upper critical field Hc2{T). Both Fermi-surface sheets were 
weighed with the same factor (1:1). The experimental result 
for LuNi2B2C is labeled with open circles. 



FIG. 3. Cp/T vs. T plot for Ya,Lui_a,Ni2B2C at zero ex- 
ternal field. 



FIG. 4. Cp/T vs. plot for Y^Lui_a;Ni2B2C at an ex- 
ternal field of B = 9 T. The straight dashed lines indicate the 
extrapolation of the data to T — > 0. 



FIG. 5. Normal-state contribution to the specific heat of 
Y^Lui_a,Ni2B2C. 



FIG. 6. Model plionon spectra obtained from the phonon 
contribution to the specific heat. The spectral densities 
a'^F{oj) used in the analysis of Ya;Lui-2!Ni2B2C were obtained 
by applying a coupling function a^{ij) ^ lu^^^^ to the model 
PDOS. The three dimensional plot shows the evolution of 
modes vs. x of F{u!). 

FIG. 7. Sommerfeld constants used in the theoretical de- 
scription compared to those obtained from specific heat mea- 
surements. 

FIG. 8. Specific heat difference AC{T) = Cs{T) - C„(T) 
(a), thermodynamic critical field Hc{T) (b), and deviation 
function D{T) (c) of LuNi2B2C compared to theoretical 
calculations for an isotropic and an anisotropic case with 
{al) = 0.03. The upper critical field (d) is a fit to the ex- 
periment with anisotropy parameters (0^) = 0.02 and 0.03. 



FIG. 9. Specific heat difference AC(T) = Cs(T) - C„(T), 
of Y2,Lui_2,Ni2B2C with x = 0.1, 0.5, 0.8, and 1 compared 
to theoretical calculations for an isotropic and an anisotropic 
case with (o^) = 0.03. 

FIG. 10. Thermodynamic critical field Hc{T) for the sam- 
ples with X = 0.1, 0.5, 0.8, and 1 compared to theoreti- 
cal calculations for an isotropic and an anisotropic case with 
(ai) = 0.03. 
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FIG. 11. Upper critical field Hc2 of Y^Lui_^Ni2B2C with 
X — 0.1, 0.5, 0.8, and 1 compared to theoretical calculations 
for anisotropic cases with (a^) = 0.02 and 0.03. 

FIG. 12. Deviation function D{T /Tc) of Y^Lui_^Ni2B2C 
with X = 0.1, 0.5, 0.8, and 1 compared to theoretical calcula- 
tions for an isotropic and an anisotropic case with {a1) = 0.03. 

FIG. 13. The critical temperature T^, electron-phonon cou- 
pling strength A, Sommerfeld constant 7, and density of states 
A^(0) as a function of x for the series Ya;Lui_a;Ni2B2C. 

FIG. 14. Specific heat difference AC(r) = C,(r) - C„(r), 
thermodynamic critical field Hc{T), and magnetic devia- 
tion function D{T) of La3Ni2B2N3_5 compared to theoret- 
ical calculations for an isotropic and an anisotropic case with 
(a^) — 0.08 . The upper critical field is a fit to the experiment 
with anisotropy parameters (a^) — 0.06 and 0.08 . 

FIG. 15. Spectral density a'^F(cj) used in the analysis 
of La3Ni2B2N3_5 obtained from a model phonon spectrum, 
which was rescaled to give the critical temperature of the sam- 
ple Tc = 11.73 K and for a fixed pseudopotential fi* — 0.13 . 

FIG. 16. Deviation functions D{T /Tc) = D{t) vs. of 
Y^Lui-_^Ni2B2C and La3Ni2B2N3_5. The BCS-resuh (solid 
line) was added for comparison. 

FIG. 17. Electronic specific heat in the superconducting 
state Cea{T) /^Tc vs. the inverse reduced temperature T^/T 
in a semilogarithmic plot. The solid line was calculated by 
adding the normal state contribution 7cT to the specific-heat 
difference between the superconducting and normal state AC 
obtaiped from the Eliashberg-equations. The data of Nohara 
et al.lB are labeled by filled triangles. 
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